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SECTION  I 


INTRODUCTION 

The  Air  Force  Flight  Dynamics  Laboratory  (AFFDL)  is  involved  in 
development  and  testing  of  techniques  to  measure  enthalpy  in  an  arc-heated 
wind-tunnel  flow.  One  of  these  methods  employs  the  spectroscopic  analysis 
of  copper  line  radiation  associated  with  an  unknown  amount  of  copper 
contaminant  that  is  present  in  the  high-temperature  effluent. 

Several  exploratory  studies  have  been  performed  using  this  non- 
obtrusive  diagnostic  technique  to  measure  stagnation  enthalpy  in  the 
AFFDL  Reentry  Nose  Tip  (RENT)  Facility^^  .  Results  of  these  efforts 
Indicated  that  a  simple  ratio  of  the  intensities  of  lines  at  5106  and  5153  A 
could  be  used  to  infer  gas  temperature  provided  that  the  following  assump¬ 
tions  are  made: 

•  Optically-thin  gas  (i.e.,  no  self-absorption) 

•  Uniform  temperature  and  copper  density  along  the 
observation  path 

•  Local  Thermodynamic  Equilibrium  (LTE) 

•  Constant  pressure. 

The  first  assumption  was  found  to  be  most  restrictive^^’ and  theoretical 
investigations  were  undertaken  at  AFFDL  to  provide  a  method  for  correcting 
the  simple  line-ratio  technique  for  self-absorption  effects^^\ 

This  report  presents  the  analytical  methods  required  to  make  the 
above  correction  and  associated  error  estimate.  In  Section  II,  the  general 
aspects  of  self-absorption  are  discussed  with  reference  to  overall  method¬ 
ology  for  accomplishing  the  correction  of  this  error  source.  Section  III 
provides  a  detailed  discussion  of  the  myriad  of  line-broadening  mechanisms 
and  associated  theoretical  models.  The  equation  forms  and  empirical 
constants  are  presented  which  provide  methods  of  calculating  the  broadening 
and  shift  of  the  copper  lines  of  interest  for  RENT  flow  conditions.  Line 
intensity  profile  theory  is  discussed  in  Section  IV  with  emphasis  on  the 
so-called  "Voight  Profile"  for  predicting  detailed  line  shape  for  varying 
copper  density.  Integrated  line-intensity  formulations  are  also  presented 
for  special  cases  of  optically-thin  and  -thick  limits.  The  computational 
methods  are  presented  in  Section  V  and  results  of  the  computer  calculations 


1 


j 

I 

I 


’1 


of  line  shape,  width,  peak  intensity,  and  integrated  line  intensity  are 
given  in  Section  VI.  Copper  density  is  treated  as  a  separate  Independent 
variable  along  with  gas  temperature,  so  the  plots  generated  in  Section  VI 
can  be  used  to  correct  measured  line-intensity  ratios  for  the  presence  of 
an  unknown  and  uncontrolled  amount  of  copper  in  the  test  gas. 

In  Section  VII,  an  error  estimate  is  provided  to  indicate  the 
accuracy  of  this  self-absorption  correction  procedure  assuming  all  other 
error  sources  are  negligible.  The  estimate  is  based  on  typical  data 
obtained  in  the  RENT  gas-cap  environment  for  peaked-profile  (l.e.,  high- 
swirl)  operation. 
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SELF-ABSORPTION  EFFECTS  ON  COPPER  SPECTRA 


For  an  optically-thin  gas  (i.e.,  no  self-absorption)  the  inte¬ 
grated  intensity  of  an  emission  line,  is  a  linear  function  of  total 

copper  number  density,  and  the  observed  length  of  the  emission  region, 

i,  such  that 


line 


"ocu  ‘  «« 


(1) 


where  T  is  the  gas  temperature.  For  this  case,  the  temperature  can  be 
determined  by  using  a  simple  ratio  of  two  emission  lines  to  eliminate  all 
other  variables  such  that 


I,  .  /I,,  =  =  f(T) 

llne^^  line2  I2 


(2) 


The  width,  y,  of  these  optically-thin  lines  is  Independent  of  copper 
density  and  optical  path  length  such  that 

Y  =  f(T,p)  (3) 


where  p  is  the  gas  pressure.  However,  as  the  copper  density  increases, 
line  radiation  emitted  by  thermally  excited  copper  atoms  from  a  particular 
point  in  the  plasma  can  be  reabsorbed  by  other  copper  atoms  in  a  lower 
energy  state  before  the  radiation  reaches  the  observer.  The  center  of  an 
emission  line  is  most  affected;  hence,  self-absorption  causes  an  increase 
in  the  measured  line  width  and  a  decrease  in  the  measured  line  peak 
intensity,  I^.  The  line  width  and  integrated  intensity  therefore  become 
nonlinear  functions  of  copper  density  and  path  length  such  that 


^llne  = 

Y  -  f(T,p,No^^,!L) 


(4) 

(5) 


For  those  transitions  associated  with  the  atomic  ground  state,  inversion 
of  the  line-intensity  profile  can  actually  occur  at  the  line  center.  But, 
for  the  particular  upper  state  transitions  of  Interest,  self-absorption 
simply  causes  a  change  in  line  contour  from  Its  original  Gaussian  shape 
to  a  flat-top  line  with  Gauss lan-shaped  line  wings. 
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Self-absorption  would  not  be  such  a  difficult  problem  if  the 
shape  and  intensity  of  all  lines  were  nearly  equally  changed.  However, 
the  more  intense  emission  lines  (i.e.,  lines  with  a  higher  strength 
factor)  are  affected  more  strongly.  For  this  reason,  the  5106  A  line 
displays  the  effects  of  self-absorption  sooner  (i.e.,  at  lower  copper 
densities)  than  the  5153  A  line.  Therefore,  the  ratio  of  the  5106  to 
5153  A  line  intensities  is  altered  by  self-absorption.  Since  this  line- 
intensity  ratio  is  used  to  infer  gas  temperature,  the  apparent  temperature 
obtained  from  the  measured  line  ratio  must  be  corrected  for  self-absorption 
effects . 

There  are  two  common  methods  of  correcting  for  self-absorption, 
and  both  involve  determination  of  the  unknown  amount  of  copper  in  the  flow 
by  measuring  another  dependent  parameter.  The  first  technique  involves 
measuring  the  absolute  intensity  of  one  of  the  lines  as  well  as  the  ratio 
of  the  two  lines  for  a  fixed  pressure  and  optical  path  length.  This  is 
an  acceptable  method  but  cumbersome  since  any  technique  involving  absolute 
intensity  measurements  is  subject  to  a  host  of  measurement  errors.  The 
method  can  be  represented  analytically  as 

p,)l  -  fixed 

=  f(T,No^^,p,Jl)  (6) 

^2^  “  f  (T»Nq^^,P,  i) 

Hence,  we  have  two  equations  and  two  unknowns  so  the  gas  temperature  and 
copper  density  can  be  determined  from  measurements  of  the  line- intensity 
ratio  12/^2  absolute  intensity  of  one  line  (i.e.,  1^^) . 

The  second  method  of  correcting  for  self-absorption  is  to  take 
advantage  of  the  fact  that  the  line  width  is  also  a  function  of  copper 
density  and  temperature  at  fixed  pressure  and  path  length.  The  set  of 
equations  can  then  be  written  whereby 


I! 


1 

i 


i 


h. 


1^' 
k ' 


I 


*  f(T.No^^.p,4)  (7) 

Yj^  “  f »No0y»P » A)  • 

Again,  the  temperature  and  copper  density  of  the  gas  can  be  determined  by 
measuring  the  llne-lntenslty  ratio  and.  In  this  case,  the  width  of  one 
line  (l.e.,  Yj^)  is  also  measured.  The  line  width,  or  full  width  of  the 
line  at  half  the  peak  Intensity,  Is  not  dependent  on  the  value  used  for 
the  absolute  magnitude  of  the  line  peak  Intensity.  Therefore,  less  error 
will  result  for  this  latter  method  of  analysis. 

In  both  methods  outlined  above.  It  Is  presumed  that  the  functional 
forms  for  the  line  width  and  line  Intensity  can  be  given  analytically  [l.e.. 
Equations  (4)  and  (5)].  These  functional  dependencies  are  required  before 
line  width  or  line  Intensity  data  can  be  used  to  correct  for  self-absorption. 
The  resulting  correction  simply  Involves  accounting  for  how  much  of  the  llne- 
lntenslty  ratio,  the  absolute  Intensity,  and  the  line  width  are  due  to 
temperature  and  pressure  effects  and  how  much  are  caused  by  the  presence 
of  excess  copper  atoms  along  the  known  optical  emission  path.  The 
analytical  representation  of  these  effects  will  be  the  subject  of  the 
following  sections  of  this  report. 


SECTION  III 


LINE  BROADENING  AND  SHIFT  THEORY 


A.  General 


Shape  and  position  of  spectral  lines  emitted  from  various  atomic 
species  in  a  plasma  have  been  a  subject  of  great  interest  to  spectroscopists. 
As  a  result,  a  good  deal  of  experimental  and  theoretical  effort  has  been 
devoted  to  the  field.  This  interest  is  motivated  by  the  fact  that  line 
width  and  shift  information,  accompanied  by  line-intensity  data,  can  be 
used  to  infer  various  plasma  properties.  In  particular,  an  accurate  measure¬ 
ment  of  stagnation  temperature  in  the  RENT  environment  cannot  be  obtained 
without  considering  line  broadening  in  order  to  correct  for  self-absorption 
effects.  However,  the  problem  of  developing  a  general  treatment  of  line 
broadening  and  shift  is  complicated  by  the  fact  that  many  mechanisms  simul¬ 
taneously  exist  which  can  cause  broadening  and  wavelength  shift  of  lines 
in  a  plasma  environment.  These  various  mechanisms  can  be  classified  as 

(1)  Fixed  broadening  and  shift 

a.  Broadening  due  to  hyperflne  structure  and 
isotope  shift 

b.  Natural  line  width 

c.  Broadening  due  to  the  measurement  device 

(2)  Doppler  broadening  and  shift 

(3)  Self-absorption  broadening 

(4)  Zeeman  broadening  in  a  magnetic  field 

(5)  Collision  broadening  and  shift 

a.  Van  der  Uaals  effect 

b.  Stark  effect  (ions  and  electrons) 

c.  Resonance  processes. 

Simultaneous  action  of  all  the  above  mechanisms  makes  it  difficult 
to  evaluate,  in  a  general  manner,  the  relative  role  of  each  process  on  the 
total  width  of  a  particular  line  of  Interest.  Further,  each  line  will  behave 
differently  even  when  the  line  originates  from  the  same  plasma  species. 

In  the  following  sections,  the  above  mechanisms  will  be  examined 
in  greater  detail  with  an  emphasis  on  the  shape  and  shift  of  the  5106  and 
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5153  A  copper  lines  In  an  arc-heated  gas  containing  an  unknown  and  uncon¬ 
trolled  amount  of  copper  contaminant.  At  first,  particular  emphasis  will 
be  placed  on  the  broadening  mechanisms  that  are  independent  of  copper 
density  in  the  flow  (i.e.,  the  optically-thin  limit).  Self-absorption 
broadening  will  then  be  the  subject  of  a  later  section  of  this  report. 


B.  Fixed  Broadening 


Fixed  broadening  is  common  to  all  spectral  lines  and  is  due  to 
natural  line  width,  hyperfine  line  structure  (i.e.,  including  isotopic 
shifts),  and  apparent  broadening  caused  by  the  measurement  apparatus.  For 
the  case  of  natural  line  broadening,  external  fields  are  neglected  (i.e.. 
Isolated  atom)  and  Heisenberg's  uncertainty  principle  is  used  to  estimate 
the  resulting  line  width  due  to  the  quantum  mechanical  uncertainty  in  the 
value  of  atomic  energy  levels.  Results  of  these  calculations  show  that 
the  natural  line  width,  Yjj»  of  the  order  of  0.0001  A.  This  width  is 
negligible  compared  with  other  line-broadening  mechanisms  and  can  be 
neglected. 


Hyperfine  structure  of  a  line  is  due  to  the  fact  that  a  given 
"line"  can  actually  be  made  up  of  several  close  but  separate  transitions. 
Nuclear  spin  effects  and  isotope  shift  cause  this  multiple  line  phenomenon; 
and  when  the  intensity  distribution  of  these  various  components  overlap, 
they  can  dominate  the  fixed  width  of  the  line.  High-resolution  measure¬ 
ments  of  the  copper  lines  of  interest  have  been  obtained  using  a  copper 

(  8) 

hollow-cathode  source'  .  The  spectra  show  that  the  5106  line  is  more 

accurately  described  as  a  triplet  of  equally-spaced  peaks,  0.02  A  apart. 

Hence,  the  5106  line  will  have  a  fixed  width  of  approximately  0.04  A.  On 

the  other  hand,  the  5153  A  line  is  truly  a  single  line  having  no  evidence 

of  hyperfine  structure  or  isotope  shift.  In  summary,  the  true  fixed  width, 

Y-,  of  the  lines  of  Interest  is  given  by 
r 


Y„(5106)  =  0.04  A 
F 

Yp(5153)  -  0.0001  A 


(8) 


7 


The  above  fixed  width  would  not  be  measured  In  practice 
because  the  apparatus  used  to  make  a  line-width  measurement  will  normally 
add  some  apparent  width  to  the  line.  This  is  due  to  the  fact  that 
the  measured  line  is  actually  a  convolution  of  the  true  line  shape  with  the 
finite  width  apparatus  function  of  the  measuring  instrument.  If  the 
Instrument  apparatus  function  has  an  inf Initely-narrow  width,  then  the 
convoluted  profile  would  represent  the  true  line  profile.  Within  certain 
limits,  apparent  line  broadening  due  to  the  experimental  apparatus  can  be 
subtracted  out  of  the  measured  line  profile  to  arrive  at  the  true  line 
width.  This  correction  procedure  involves  determining  the  apparatus  func¬ 
tion  experimentally,  then  deconvolutlng  the  measured  profile  with  the 
apparatus-function  profile  to  obtain  the  original  line  contour.  The 
apparatus  function  is  measured  by  using  that  instrument  to  scan  a  light 
source  that  emits  radiation  with  ultra-narrow  lines  (i.e.,  source  lines 
are  approximately  100  times  narrower  than  apparatus-function  width) . 

Lasers,  pen-ray  lamps,  or  hollow-cathode  light  sources  produce  such  ultra¬ 
narrow  lines. 

Provided  that  the  apparatus  function  and  the  measured  line  profile 
have  characteristic  shapes  (i.e.,  Gaussian,  triangular,  etc.),  the  convolu¬ 
tion  can  be  performed  once  for  various  full  half-widths  and  then  correction 
curves  can  be  used  from  that  point  on  to  infer  the  true  line  width  from  the 
measured  width  and  the  apparatus-function  width.  Such  correction  curves 
for  several  characteristic  line  shapes  are  shown  in  Figure  1.  For  the 
particular  conditions  of  the  present  experimental  setup,  both  the  true  line 
and  apparatus-function  profiles  are  nearly  Gaussian  shaped.  The  true  line 

width,  Y,  can  then  be  determined  from  the  measured  width,  y  ,  by  using  the 

m 

solid  curve  in  Figure  1  given  a  measured  value  of  the  apparatus-function 
width  (i.e.,  y^^) •  Note  that  the  correction  for  apparent  line  broadening 
becomes  quite  large  for  Y^/y^  >  0.5.  Therefore,  an  effort  should  be  made 
to  keep  the  apparatus  function  as  narrow  as  possible  to  avoid  errors  in 
correcting  for  apparent  line  broadening.  Below  the  value  of  y.  =  y  /3,  the 
measured  width  is  very  close  to  the  true  width  (i.e.,  y  ==  y) . 


FIGURE  1.  CORRECTION  FACTOR  TO  ACCOUNT  FOR  APPARENT  LINE 
BROADENING  DUE  TO  THE  MEASUREMENT  APPARATUS 


C.  Doppler  Broadening 

The  steady  motion  of  a  radiating  particle  along  the  direction 
of  sight  of  an  observer  leads  to  a  discrete  wavelength  shift  of  the  line 
which  is  called  a  Doppler  shift.  Furthermore,  in  a  high-temperature  plasma, 
the  random  motion  of  the  particles  toward  and  away  from  the  observer 
produces  a  smearing  out  of  the  emitted  line.  This  motion- induced,  line- 
broadening  phenomenon  gives  rise  to  the  Doppler  line  width.  In  our  case, 
the  high-speed  arc-heated  flow  is  being  observed  at  an  angle  perpendicular 
to  the  gas  flow  direction.  Hence,  only  Doppler  broadening  due  to  the  random 
motion  of  the  radiating  particles  occurs  (i.e.,  no  discrete  flow-induced 
Doppler  shift  is  observed). 

If  the  component  of  velocity  along  the  line  of  sight  of  the 

observer  is  v  ,  then  the  discrete  shift  of  the  line,  L\  ,  due  to  the  Doppler 
s  s 

effect  can  be  written  as 

V 

AAg  =  +  -^  X  [A]  (9) 

where  A  is  the  unshifted  wavelength  of  the  line  and  c  is  the  velocity  of 
light.  Assuming  that  the  plasma  particles  move  in  a  random  manner  due  to 
thermal  energy,  the  distribution  of  velocities  are  Maxwellian  in  nature. 
Hence,  the  fraction  of  radiating  particles,  dN/N,  in  the  observer's  line 
of  sight  with  velocity  components  between  v^  and  v^  +  Av^  is 

—  =  —  exp  t“(v  /v  )  ]  -:r—  (10) 

—  sp  v 

V  Tl  ^  P 

where  v^,  the  most  probable  velocity  for  a  particle,  is  given  by 


In  Equation  (11),  R  is  the  universal  gas  constant,  T  is  the  gas  temperature, 
and  is  the  atomic  weight  of  the  emitting  atom.  Substituting  Equation  (9) 
into  Equation  (10)  yields 

- - -  exp  [-(AA/AA  )^]  d(AA)  (12) 

N  A"  AA  P 

P 
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where  now  AA  is  defined  by 
P 


(13) 


For  optically-thin  conditions,  the  line  intensity  will  be  proportional 
to  the  concentration  of  radiating  particles,  N.  Hence,  the  Intensity 
emitted  in  the  interval  I(AA)  d(AA)  is  proportional  to  the  fraction  of 
particles  dN  given  in  Equation  (12).  Therefore,  the  line  intensity  for  a 
Doppler-broadened  line  can  be  written  as 


I(AA)  =  - ^ —  exp  [-(AA/AA  )^] 

/n  AA  ^ 

P 


(14) 


where  is  the  total  integrated  line  intensity.  The  form  of  Equation  (14) 
indicates  that  broadening  due  purely  to  the  Doppler  effect  results  in  a 
Gaussian  line  shape.  The  full  width  of  this  line  (i.e.,  the  total  width 
between  points  where  the  intensity  falls  to  half  its  value  at  the  maximum) 
is  given  by 


1/2 

=  2(ln2)  '  AA  =  1.66  AA 

D  p  p 


(15) 

Yp  =  7.16  X  10"^ 

A  in  A,  T  in  "K,  y  in  AMU 

where  Yq  is  defined  as  the  Doppler  full  width  of  the  line.  From 
Equation  (15) ,  it  can  be  seen  that  this  effect  is  most  pronounced  for 
light  atoms  at  high  temperature.  For  copper,  y^  =  63.56,  Equation  (15) 
yields 

Yjj(copper)  =  8.98  x  lO"®  A  [A]  .  (16) 

The  line  width  for  the  copper  lines  near  5100  A,  and  a  temperature  of 
Interest  (i.e.,  6000  K) ,  is  given  by 

Yp(5106)  =<  Yo(5153)  -  0.036  A  .  (17) 

This  width  is  of  the  same  order  as  the  fixed  broadening  value  for  the 
5106  line. 
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D.  Zeeman  Broadening 


In  the  presence  of  a  magnetic  field,  an  atomic  line  splits  into 
several  lines.  This  splitting  is  small  and  symmetric  about  the  original 
line  position.  The  resulting  effect  is  to  produce  a  broadening  of  the 
original  line  when  the  multiplet  is  viewed  under  low  resolution.  Exact 
calculations  are  complex,  but  an  estimate  can  be  made  which  is  fairly 
accurate.  This  estimate  indicates  that  the  line  full-width,  Zeeman 

broadening  in  a  magnetic  field  of  strength  B  can  be  written  as 

=  10"^  B  [A]  (18) 

X  in  A,  B  in  kllogauss. 

In  the  RENT  facility,  two  coils  are  used  to  rotate  the  arc 
attachment  point  in  order  to  extend  electrode  life.  The  coil  nearest  the 
measurement  point  has  a  magnetic  field  strength  of  2  kilogauss  at  its 
center.  However,  at  the  measurement  point,  the  field  strength  has  decreased 
to  approximately  0.1  kilogauss.  Hence, 

Y  (5106)  =  Y  (5153)  =  0.0026  A  .  (19) 

z  z 

This  is  a  small  effect  and  can,  therefore,  be  neglected  for  the  present 
RENT  conditions. 


E.  Collision  Broadening  and  Shift 

In  this  case,  broadening  and  shift  of  a  line  occur  due  to  the 
interaction  of  radiating  atoms  with  surrounding  particles.  For  a  plasma 
environment,  these  Interacting  particles  could  be  electrons  or  ions  (i.e.. 
Stark  broadening),  neutral  molecules  or  nonsimilar  atoms  (i.e..  Van  der  Waals 
broadening),  and  atoms  of  the  same  kind  as  the  emitting  atom  (i.e.,  resonance 
broadening).  In  the  present  case,  radiation  from  copper  contaminant  in 
arc-heated  air  is  being  considered  so  copper-copper  Interactions  are  rare. 
Hence,  resonance  broadening  will  be  neglected. 

The  analytical  representation  of  collision  broadening  is  the 
most  difficult  of  all  the  broadening  mechanisms  because  it  involves  the 
modeling  of  atom-particle  interactions.  Furthermore,  the  overall  broadening 


of  a  particular  line  depends  on  the  type  of  atom  and  the  particular  energy 
states  Involved  in  the  radiating  transition.  The  formulation  also  depends 
on  the  particular  type  of  Interacting  partner.  Hence,  much  more  experi¬ 
mental  e'v’ldence  is  required  to  obtain  empirical  broadening  constants  and 
to  validate  the  theory  for  a  particular  line  of  interest. 

The  analytical  solution  to  the  problem  of  predicting  collision 
broadening  can  be  set  up  on  the  basis  of  the  quantum-mechanical  theory  of 
atomic  collision.  However,  the  results  of  quantum  theory  calculations 


nearly  always  agree  with  simple  formulations  that  can  be  obtained  using  a 

.(9) 


quasl-classlcal  approach  .  This  latter  approximation  considerably 
simplifies  the  closed  form  solution  to  the  problem  and  makes  the  process 
of  evaluation  more  physically  meaningful. 

In  quasi-classical  collision  theory,  the  interaction  of  the 
emitting  atom  with  surrounding  particles  is  considered  as  an  amplitude- 
phase  modulation  of  the  unperturbed  oscillator.  The  Intensity  distribution 
in  an  emission  line  is  then  determined  by  performing  a  Fourier  analysis 
averaged  over  all  collisions.  However,  in  performing  the  mathematical 
computations,  certain  inequalities  must  be  fulfilled  and  this  leads  to 
several  different  regimes  of  interest.  It  can  be  shown  that  the  closest 


particle-atom  interaction  creates  the  greatest  perturbation  and  therefore 

.(10) 


contributes  mostly  to  the  wings  of  the  line  contour  '.  The  combined 
action  of  the  rest  of  the  particles  affects  the  line  center.  The  closest 
particle  interactions  are  considered  best  represented  by  a  general 
statistical  theory  approach  while  the  interactions  due  to  longer  range 
encounters  are  amenable  to  the  so-called  impact  theory  methods.  It  should 
be  noted  that  neither  theory  can  be  strictly  applied  over  the  whole  line 
profile.  In  the  wings,  the  statistical  theory  is  always  valid;  and  at  the 
line  center,  the  impact  theory  holds  true.  The  impact  theory  computations 


are  simplified  by  assuming  that  the  atom-particle  interaction  time,  Tj,  is 


short  compared  to  the  time  between  collisions,  (i.e.,  Xj  <<  Tj,)  •  The 


impact  theory  is  therefore  strictly  valid  only  for  low  pressures.  On  the 
other  hand,  the  statistical  theory  is  strictly  valid  only  for  very  high 
pressures. 


13 


:  i. 


■.  1 


V 


I 


1.  Impact  Theory 


As  discussed  above,  the  intensity  distribution  near  the  line 
center  is  represented  by  the  so-called  impact  theory.  The  basic  principles 
behind  this  theory  of  line  broadening  were  proposed  by  Lorentz.  He  assumed 
the  emission  of  an  atomic  oscillator  is  interrupted  during  a  collision  with 
perturbing  particles.  After  a  collision,  emission  at  the  same  frequency 
resumed  with  no  emission  during  the  interaction.  Hence,  vibration  of  the 
oscillator  is  represented  by  a  series  of  sinusoidal  segments  of  duration  T 
(see  Figure  2a) .  A  Fourier  analysis  of  such  oscillations  gives  the  following 
expression  for  the  line  intensity: 


Ko)) 


_1 _ 1 _ 

’'^r  (o)^  -  co)^  +  (1/t^)^ 


(20) 


where 


10  =  27rv  .  (21) 

In  the  above  equations,  v  is  the  frequency  of  radiation,  the  mean  radi¬ 
ation  lifetime,  and  (o^  the  circular  frequency  of  the  line  center.  The 
intensity  is  normalized  such  that 


CD 

/  1(a))  d  0)  =  1  .  (22) 

00 


Equation  (20)  is  the  so-called  Lorentz  dispersion  formula  for  collision 
broadening.  The  shape  of  the  line  is  often  called  a  Lorentz  shape  or 
dispersion  shape.  The  collision-broadening  constant,  y is  defined  by 


Y 


c 


(23) 


which  represents  the  full  width  of  the  line  at  half  the  maximum  intensity 
in  units  of  circular  frequency.  Now,  can  be  related  to  the  collision 
frequency,  z,  if  we  assume  that  the  radiative  lifetime,  x^,  is  of  the  order 
of  the  time  between  collisions,  x^.  Hence, 


z 


(24) 
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(a) 


FIGURE  2,  VIBRATION  OF  A  PERTURBED  OSCILLATOR 
(a)  LORENTZ  (b)  LENZ-WEISSKOPF 
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Thus,  by  defining  an  effective  collision  cross  section,  Tfp  „  ,  the  colli- 
sional  half-width,  y^>  can 
collision  theory  to  obtain 


sional  half-width,  y^>  can  be  related  to  plasma  properties  using  ordinary 


Y  =  ~  =  2Trp  ^  N  V  [sec  (25) 

CT  CPD*"  ' 

C  ft' 

where  is  the  number  density  of  the  perturbing  particles,  is  the 
effective  collision  radius,  and  v^,  the  relative  collision  velocity,  is 
given  by 


.  /MI  (i  ^  J_, 


(26) 


In  Equation  (26),  and  ^2  ^he  weights  of  the  colliding  partners  in 
atomic  mass  units. 

Since  the  number  density  of  perturbers  is  proportional  to  the 
gas  pressure,  the  line  broadening  given  in  Equation  (25)  would  be  propor¬ 
tional  to  pressure  and,  hence,  the  name  "pressure  broadening"  is  sometimes 
used, 

A  shortcoming  of  the  Lorentz  formulation  is  that  it  does  not 
define  a  means  of  computing  the  effective  collision  radius,  p^,  for  various 
types  of  interactions  (i.e.,  atom-electron,  atom-neutral,  etc.).  It  also 
does  not  account  for  line  shift  due  to  the  atom-particle  interaction. 

To  rectify  the  shortcomings  of  Lorentz's  simple  model,  Lenz' 

postulated  that  the  oscillator  did  not  stop  emitting  during  a  collision, 

but  rather,  a  change  in  frequency  of  oscillation  occurred  during  a  collision 

(see  Figure  2b) .  Hence,  the  colliding  particles  could  also  contribute  to 

the  broadening  of  a  line  by  a  so-called  "phase  altering"  collision  without 

extinguishing  the  radiation.  In  this  interaction  process,  as  visualized 

(12) 

by  Lenz  and  Weisskopf  ,a  collision  acts  to  detune  (i.e.,  broaden  the 
frequency  spectrum)  of  a  radiator. 

Following  this  model,  a  particle  at  a  distance  r  from  a  radiating 
atom  causes  a  phase  change.  Am,  of  the  form 


(11) 


2Trc 


Am  = 


(27) 
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where  the  exponent  n  depends  on  the  particular  type  of  force  law  which 

represents  the  interaction  between  the  perturbing  particle  and  the  radiating 

atom.  In  the  case  of  the  Stark  broadening  effect,  this  exponent  has  the 

value  of  4  and  for  a  Van  der  Waals  interaction,  n  =  6  is  used.  The  constant 

c  characterizes  the  interaction  of  the  electric  field  of  the  perturbing 
n 

particle  with  the  particular  energy  level  of  the  radiating  atom.  Hence, 
c^  can  vary  widely  for  each  atomic  emission  line  and  for  each  type  of 
perturbing  particle. 

In  this  formulation,  the  distance  of  closest  approach  to  the 
radiating  atom  is  designated  by 

Vn  =  ^  •  ^28) 

Results  of  the  computations,  following  Lenz  and  Weisskopf, 
show  that  the  total  phase  change,  n,  due  to  collision  with  a  particle 
having  an  impact  parameter  p^,  is  given  by 


solving  Equation  (29)  in  terms  of  the  impact  parameter  we  ge:. 


j-n^  c^\  l/(n-l) 


According  to  Lenz  and  Weisskopf,  an  effective  collision  (i.e.,  one  that 
causes  line  broadening)  occurs  for  >  n  “  !•  Therefore,  divides  all 
collisions  into  two  kinds 


n  <  tIq  -  no  broadening 


n  >  -  broadening 


Lenz  and  Weisskopf  arbitrarily  picked  n 
collision  radius  of  the  form 


1,  which  gives  an  effective 


l/(n-l) 
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where  Is  often  called  the  Uelsskopf  collision  radius.  Also,  Is 
sometimes  referred  to  as  the  optical  radius.  The  full  half-width  In 
circular  frequency  units  for  the  Lenz-Weisskopf  collision  theory  is 
therefore  given  by 


Y 


c 


V 

P 


(33) 


where 


Stark  broadening 


Van  der  Waals  broadening 


(34) 


The  limit  of  n  =  =  1  in  Lenz's  theory  is  an  approximation 

which  results  in  the  fact  that  only  certain  collisions  are  taken  into 

account.  This  approximation  limits  Lenz's  theory  in  principle  to  very 

low  pressures  and  to  the  intensity  distribution  at  the  core  of  a  line. 

However,  the  theory  does  predict  a  possible  line  shift  as  well  as  line 

broadening  and  provides  a  method  to  compute  effective  collision  radii  for 

various  t3rpes  of  atom-particle  interactions. 

(13) 

Lindholm  succeeded  in  overcoming  the  difficulties  in  the 
Lenz-Weisskopf  theory  by  placing  no  limit  on  the  magnitude  of  the  phase 
shift  and  by  taking  all  phase  shifts  into  account  (i.e.,  both  close  and 
distant  encounters  are  considered).  Lindholm' s  more  general  theory  is 
valid  at  higher  pressures  and  over  a  greater  portion  of  the  line  profile. 
It  also  retains  the  good  points  of  Lenz's  approach.  Detailed  calculations 
show  that  the  intensity  distribution  is  of  the  form 


Km) 


X. 

2ir 


_ 1 _ 

(m  -  -  A)^  +  (y/2)^ 


(35) 


where  the  intensity  has  again  been  normalized  to  unity  as  in  Equation  (22). 

In  Equation  (35),  A  is  the  line  shift  in  circular  frequency  units.  Lindholm' 
analysis  resulted  in  the  similar  formulation  for  the  effective  collision 
radius  and  line  full  half-width  as  given  in  Equations  (30)  and  (33),  respec¬ 
tively.  However,  details  of  the  calculations  show  that  the  collision  radius 
is  actually  given  by 
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l/(n-l) 


(36) 


a  C  \ 
n  n  \ 

*^0  *  V  n  1 

\P  oj 

where  in  Equation  (36)  is  usually  set  equal  to  0.64  for  Lindholm's 
theory  rather  than  the  arbitrary  value  of  =  1  for  Lenz's  formulation. 
The  values  of  are  the  results  of  detailed  Fourier  analyses  of  the  col¬ 
lision  model  and  empirical  constants  determined  by  appropriate 

line  width  measurements  at  known  conditions. 

Results  of  Lindholm's  calculations  also  showed  that  the  line 
shift.  A,  could  be  written  as 

A  “  p  ^  N  V  .  (37) 

o  p  p 

Hence,  from  Equations  (33)  and  (37),  Lindholm's  approach  implies 
a  constant  magnitude  of  the  ratio  of  width  to  shift  which  characterizes  a 
definite  type  of  interaction.  Detailed  calculations  show  that  according 
to  Lindholm's  theory  the  width-to-shif t  ratio  is  given  by 


=  1.16,  (Stark) 
4 


—  =  2.8,  (Van  der  Waals) 

6 

where  the  subscripts  4  and  6  refer  to  Stark  and  Van  der  Waals  broadening, 
respectively. 


2.  Validity  of  Impact  Theory  for 
RENT  Environment 


In  this  section,  the  applicability  of  utilizing  impact  theory 
rather  than  statistical  theory  to  predict  line  broadening  is  discussed  for 


the  RENT  environment.  Sobel'man' 


has  shown  that  the  collision  mechanism 


is  valid  for  predicting  the  intensity  distribution  near  the  line  center 
provided  that  a  certain  criteria  on  the  phase  shift.  Am,  is  met  such  that 

n/n-1 

^n 


t 


For  the  present.  Equation  (39)  can  be  simplified  using  Equations  (33)  and 
(36)  to  write  the  inequality  in  terms  of  the  number  density  of  perturbing 
particles  and  the  optical  collision  radius.  The  resulting  criterion  is 


(40) 


The  Weisskopf  radius  for  the  copper  lines  of  interest  is  of  the  order  of 
5  Hence,  from  Equation  (40) 


N  «  10^^  cm' 
P 


(41) 


must  be  satisfied  for  the  Impact  theory  to  be  valid  in  the  case  of  collision 
broadening  of  the  copper  lines  of  interest.  Calculations  show  that  for 
typical  RENT  conditions  (i.e.,  =  100  atm,  T^  =  6000  K)  number  density 

levels  in  the  free  stream  and  behind  the  bow  shock  wave  are 


Np(RENT) 


22  -3 

«  10  cm 


Hence,  typical  RENT  flow  conditions  are  such  that  the  impact  theory  can 
be  considered  valid  for  predicting  line  broadening  near  the  line  center. 
Further  evidence  of  the  applicability  of  impact  theory  will  also  be  provided 
when  measured  line  shape,  full  half-width,  and  shift  are  compared  to 
theoretical  predictions.  Details  of  the  impact  theory  are  presented  in 
the  following  text. 


3.  Modifications  to  Lindholm* s 


Impact  Theory 


The  applicability  of  the  quasl-classical  theory  of  Lindholm  to 
electron-atom  interactions  (i.e..  Stark  broadening)  is  more  complex  than 
the  above  analysis  indicates.  This  is  due  to  the  fact  that  for  electrons 
the  Weisskopf  mechanism  for  line  broadening  becomes  less  effective  due  to 
the  presence  of  inelastic  collisions.  Interactions  of  this  type  (i.e., 
those  that  involve  energy  exchange  between  the  collision  partners)  were 
neglected  in  Lindholm' s  semi-classical  theory.  More  complete  analytical 
results  have  been  performed  by  Vainshtein  and  Sobel'man^^^^  which  Include 
these  inelastic  effects.  The  resulting  formulation  Indicates  that 
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(43) 


f«4  S\l/3 


n  V 
o  e 


l4  ,  ,  16  l'(3) 

1.16  j„(3j 


where  I' (3).  I"(3)  are  correction  factors  that  are  dependent  on  the 
elasticity  parameter  3  which  varies  for  each  line  of  interest.  The  3 
parameter  can  be  written  as  follows: 


3 


<r 


3/2 


(44) 


where  AE  is  the  distance  to  the  nearest  energy  level  that  can  exchange 
energy  with  the  original  level  involved  in  the  emission  line.  This  energy 
exchange  between  levels,  which  is  facilitated  by  low  AE's,  leads  to  a 
higher  degree  of  inelastic  interaction  during  atom-electron  collisions. 

For  lines  with  large  3  values  (i.e.,  large  AE's  or  3  >  1) ,  I' 
and  I"  are  approximately  1,  so  no  correction  results.  However,  for  small 
3's,  I'  and  I"  are  less  than  1  with  I"  falling  off  faster  than  I'.  Hence, 
the  inclusion  of  inelastic  effects  will  result  in  slightly  narrower  lines 
and  a  wldth-to-shift  ratio  greater  than  1.16  for  the  Stark  broadening  case. 

Van  der  Waals  broadening  and  shift  are  not  affected  by  the  above 
correction  factors  and  can  be  written  as  follows: 


^6  =  Njj  Vj, 


“6  ^6  ' 


n  V 
o  n 


(45) 
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4.  Ion  Broadening 


In  a  plasma,  the  line  broadening  due  to  the  Stark  effect  (l.e., 
Y^)  Is  determined  by  the  Interaction  of  an  emitting  atom  with  ions  as  well 
as  electrons.  For  ions,  the  Inelastic  collision  correction  is  absent  for 
cases  of  practical  Interest  because  the  ion  velocity  is  small  compared  to 
that  of  an  electron.  In  general,  we  can  write 

^4  “  Y^(electrons)  +  y^dons)  (46) 


where  from  impact  theory  [l.e..  Equation  (43)] 

1/3 

Y^(lons)  «  Vj.  Nj 


(47) 


where  Nj  and  Vj  are  the  ion  number  density  and  velocity,  respectively. 
For  a  neutral  plasma,  =  N^,  and  the  ion  velocity  can  be  expressed  in 
terms  of  the  electron  velocity  using  Equation  (26)  to  write 


V 


I 


(48) 


where  and  Pj  are  the  electron  and  ion  mass,  respectively.  Using  this 
expression,  we  get  the  following: 


1/3  „  1/3  „ 

Vt  N,  “  V  N 
I  I  e  e 


Ue  U/6 


Hence,  using  Equation  (43),  Equation  (46)  can  be  written  as 

’v  \l/6 

Y.  “  v^  N^[I'(e)  + 


(49) 


(50) 


where  the  correction  for  ion  broadening  is  of  the  form 


v"'/ 


1/6 


(51) 


More  exact  calculations  given  in  Reference  8  show  that  the  actual  correc¬ 
tion  factor  for  ion  broadening  is  of  the  form 


(52) 
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where  Z  Is  the  ionic  charge.  Hence,  the  Stark  broadening  of  a  line 
including  inelastic  electron  collisions  and  ion  effects  is  given  by 


27TP  V  N  [I'(6)  + 
o  e  e 


2m  Z 
e 


2  \  1/6 


] 


'^e  % 


1/3 


and 


1.16 


I'(3)  + 

f  2Pg  1 

|l/6 

f2p  Z^  ^ 

I"(6)  + 

1  ^ 

i  i 

^1/6 

(53) 


For  RENT  conditions,  thermodynamic  calculations^^ show  that  for 
p  =  80  atm  and  T  =  6000  K  (i.e.,  model  gas-cap  conditions)  the  dominant 
ion  is  NO^  since 


Therefore, 


NO"* 


=  N 


z  =  1 


(54) 


Pi  “  PyQ+  =  30.02  AMU  (55) 

p  =  5.49  X  10  ^  AMU 
e 


and 


2^  \  1'* 


0.18 


(56) 


Examination  of  Equations  (53)  and  (56)  indicates  that,  for  the 
RENT  environment,  NO  ion  collisions  will  Increase  the  total  Stark  broaden¬ 
ing  by  at  least  20  percent.  However,  the  Stark  shift  will  be  affected  by 
ion  collisions  only  for  low  6  lines  (i.e.,  I"  <<  I')  since  the  term  of 
Equation  (56)  appears  in  both  the  numerator  and  demoninator  of  Equation  (53). 
Other  ions  will  not  contribute  significantly  to  line  broadening  in  the  RENT 


/ 
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environment  since  their  number  densities  are  at  least  two  orders  of 
magnitude  less  than  that  of  the  NO^  ion.  With  the  above  modification, 
Llndholm's  formulation  was  adopted  for  the  present  Impact  theory  analysis. 
Other,  more  complex  collision  theories  could  be  used.  However,  Llndholm's 
representation  does  predict  the  main  features  of  atomic  line  broadening 
and  shift  without  undue  mathematical  complication.  The  remaining  problem 
Is  to  specify  the  constants  necessary  to  calculate  line  broadening  and 
shift  for  the  RENT  flow  environment  using  Llndholm's  analytical  formulation. 

5.  Determination  of  Collision-Broadening 
Constants 


There  are  specific  constants  given  in  Equations  (45)  and  (53) 

[i.e.,  a  ,  C  ,  6,  I' (6),  I"(B)]  which  must  be  determined  for  each  line  of 
n  n 

Interest  before  collision  broadening  can  be  accurately  calculated.  Some 
of  these  constants  are  available  from  the  literature  and  others  must  be 
estimated.  It  is  the  purpose  of  this  section  to  Identify  the  best  available 
values  for  these  constants  in  order  to  compute  line  broadening  in  the  RENT 
environment. 

Table  1  gives  the  "best  estimate"  for  values  of  a^,  n^,  and 

corresponding  that  are  reported  by  various  authors.  The  term  "best 

estimate"  is  used  here  because  it  is  often  difficult  to  determine  exactly 

what  an  author  used  due  to  differences  in  defining  the  line  width  (i.e., 

full  half-width  versus  half  half-widths)  and  also  because  of  differences 

in  units  (i.e.,  2tt  factor  involved  in  converting  circular  frequency  to 

ordinary  frequency).  These  differences  are  not  really  that  important 

since  changes  in  the  empirical  constant,  C^,  which  is  experimentally 

determined,  will  correct  for  differences  in  an  author's  particular  choice 

or  definition  of  a  and  n  . 

n  o 

The  constants  for  copper  lines  of  interest  are  determined  from 
measurements  of  line  width  and  shift  at  known  thermodynamic  conditions. 
However,  a  problem  arises  in  the  fact  that  little  data  actually  exist  for 
copper  line  broadening  and  shift.  Furthermore,  most  of  the  data  that  does 
exist  comes  from  measurements  obtained  in  low  pressures  (i.e.,  1  atm)  free- 
burning  arcs^^^^  where  thermod3rnamlc  conditions  are  often  not  well 
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quantified.  AFFDL  has  sponsored  experimental  efforts  to  measure  the 

broadening  and  shift  constants  of  copper  lines  In  air  using  a  shock  tube 

operating  at  conditions  typical  of  the  RENT  flow  environment.  However, 

only  Van  der  Waals  data  (l.e.,  C,)  for  the  5106  A  line  have  been  successfully 
(18)  ** 

obtained  to  date'  .  Quantum-mechanical  methods  have  also  been  used  to 
(19) 

calculate  .  Values  of  the  collision-broadening  constants  from 

various  sources  are  given  in  Table  2.  Note  that  the  sign  of  the  C,  and  C, 

4  o 

constants  Indicates  the  direction  of  shift  of  the  line;  i.e.,  a  negative 
sign  implies  a  red  shift.  Examination  of  Table  2  indicates  that  the  majority 
of  information  on  copper  line  broadening  and  shift  comes  from  Ovechkin's  data 
from  a  free-burning  arc.  Mlyachi's  calculations  of  for  the  5153  A  line 
disagree  with  Ovechkin's  measurements  by  a  factor  of  3.  Calspan's 
results  at  elevated  pressures  (i.e.,  90  atm)  seem  to  disagree  with  Ovechkin's 
Cg  results  but  only  by  a  factor  of  1.7.  However,  Ovechkin's  formulation 
for  Van  der  Waals  broadening  uses  a  different  value  than  that  used  to 
determine  using  Calspan's  data  since  the  latter  used  the  formulation 
of  Miyachi  (see  Table  1).  Hence,  the  actual  Calspan  combination 

predicts  a  factor  of  2.6  greater  line  width  for  Van  der  Waals  broadening 
than  was  inferred  from  Ovechkin's  results.  A  decision  is  required  as  to 
which  values  given  in  Table  2  are  to  be  used  in  the  present  analysis. 

It  was  decided  that  Calspan's  Cg  value  for  the  5106  A  line  was 
more  representative  of  broadening  conditions  in  the  RENT  environment.  The 

O 

value  given  by  Ovechkin  for  C^  of  the  5153  A  line  was  then  modified  to 
reflect  the  same  order-of-magnitude  difference  in  as  measured  by  Ovechkin, 
so  that 


Cg(5106)=  -1.7  X  10 
Cg(5153)=  -  1.7  X  10 


(57) 


Since  no  other  data  are  presently  available,  Ovechkin's  values  were  assumed 
for  the  remainder  of  the  constants  given  in  Table  2. 

Up  to  this  point,  the  formulations  for  collisional  line  width 
and  shift  have  been  given  in  circular  frequency  units.  To  convert  these 
values  to  wavelength  units  in  Angstroms,  we  must  use  the  conversion  factor 
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BROADENING  CONSTANTS  REPORTED  BY  VARIOUS  AUTHORS  FOR  COPPER  LINES  OF  INTEREST 


5153 


Y.A  [A]  =  y.A  [sec  ]  X  x  10 

where  X  is  in  cm,  and  c,  the  velocity  of  light,  is  in  cm/sec. 

In  summary,  we  can  write  the  entire  set  of  equations  necessary 
to  compute  collisional  line  broadening  using  Tables  1  and  2,  and  Equa¬ 
tions  (45),  (53),  and  (57)  to  obtain 

Y,  -  17.0  V  3/5  Nn  ic  ^  ““  iSl 

6  6  N 


IT  =  2.8 


/Srt  ,  1  ^  1  ^ 


-1.7  X  10~^^,  A  =  5105.5  X  10  ® 
-1.7  X  10~^^,  X  =  5153.2  X  lO"® 


[cm  /sec] 


\-n 


■4C,3/5  v^l/3  [1.(6)+  U 


Y,  I '(3)  +  (2U  Z  /M  ) 

^  =+  1.16 - ^ ^ 

A,  -  -  J 

I" (3)  +  (2y^z^/h  ) 

e  1 


where 


I’ (3) 


I"  (3) 


5105.5  X  10 
5153.2  X  10‘ 


(cm^/ sec) 


-1.5  X  10 
8.1  X  lO' 


R||||PI4  . . .  . ■  .•  ..11.  I..  ,. . .  IIBU  I.IIPIM  »■ . 

The  above  equations  can  be  reduced  to  simpler  forms  by  substi¬ 
tuting  appropriate  values  for  the  various  parameters  to  obtain 

r,  -  5.96  X  10-“  «  I^/lO  111  (5106) 

O  N 

Y,  =  1.52  X  10~^^  N„  [X]  (5153) 

o  N 

2.8 

9.66x10-1%  ^1/6  111  (5106) 

e 

3.29  X  10“^®  N  (5153) 

e 

^  =  -1.16  (5106)  (59) 


^  =  -H2.28 


(5153) 


^c  = 


Ya  +  Y/ 


A 

c 


^6  \ 


or 


A 

c 


2.8  1.16 


(5106) 


A 

c 


2.8  2.28 


(5153) 


Results  of  calculations  using  these  equations  to  compute  line 
widths  and  shifts  due  to  collision  effects  are  shown  in  Figures  3  and  4. 

A  typical  RENT  nose-tip  pressure  of  80  atm  was  used  in  the  analysis  to 
infer  the  number  density  of  neutrals,  N^j,  and  electrons,  N^,  from  equilibrium 
thermodynamic  calculations.  Note  that  for  conditions  characteristic  of  the 
RENT  nose-tip  environment  (i.e.,  p  =  80  atm,  T  =  6000  K) ,  Van  der  Waals 
broadening  and  shift  dominate  the  Stark  effect  for  the  copper  lines  of 
Interest . 
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FIGURE  3.  COLLISIONAL  FULL  HALF-WIDTH  OF  APPROPRIATE 

COPPER  LINES  AS  A  FUNCTION  OF  GAS  TEMPERATURE 
(p  “  80  atm) 
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FIGURE  4.  LINE  SHIFT  DUE  TO  COLLISIONS  FOR  APPROPRIATE 
COPPER  LINES  AS  A  FUNCTION  OF  GAS 
TEMPERATURE  (  p  -  80  atm) 
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A, 


F.  Summary  of  Broadening  and  Shift 


In  order  to  identify  the  important  broadening  mechanisms  in  the 
RENT  environment,  a  calculation  was  made  of  the  magnitude  of  various 
contributions  to  5106  and  5153  A  line  widths  for  conditions  typical  of  an 
arc-tunnel  test.  The  tunnel  stagnation  conditions  were  assumed  to  be  as 
follows: 

p  “  100  atm 

*^o 

T  =  6000  K 
o 

From  thermodynamic  and  fluid  dynamic  calculations,  the  following  independent 
variables  were  used  in  the  analysis  for  line  broadening  behind  the  bow  shock 
wave: 

p  =  80  atm 
T  =  6000  K 

N„  =  N  =  9.66  X  10^^ 

N  o 

N  =  6.6  X  10^^ 
e 

Utilizing  formulations  of  the  previous  sections,  these  independent  vari¬ 
ables  result  in  the  computation  of  broadening  and  shift  values  found  in 
Table  3. 

Note  that  the  Van  der  Waals  effect  (i.e.,  atom-neutral  collisions) 
was  the  dominant  broadening  and  shift  mechanism  for  the  RENT  environment 
assximlng  optically-thin  lines.  This  is  fortunate  since  the  more  accurate 
Calspan  shock-tube  data  implied  Van  der  Waals  constants  only.  From  the 
results  in  Table  3,  it  was  determined  that  only  Doppler,  Van  der  Waals,  and 
Stark  mechanisms  must  be  considered  for  line- intensity- ratio  formulations 
applicable  in  the  RENT  environment. 

If  self-absorption  is  important,  the  measured  line  widths, 
corrected  for  apparatus  function  broadening,  will  be  wider  than  indicated 
in  Table  3.  Therefore,  the  above  broadening  calculations  will  be  required 
to  determine  the  unknown  copper  density  from  the  measured  line-width  values. 
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SUMMARY  OF  BROADENING  AND  SHIFT  VALUES  FOR  RENT  CONDITIONS 


Resonance 


SECTION  IV 

LINE  INTENSITY  THEORY 


A.  General  Analytical  Form 


Use  of  the  simple  llne-intenslty-ratio  technique  Implicitly 
assumes  chat  the  gas  under  observation  emits  but  does  not  absorb  radiation 
at  the  wavelengths  of  interest  (i.e.,  the  gas  is  optically  thin).  For  the 
case  where  this  assumption  Is  not  valid,  a  more  complex  formulation  must 
be  developed  to  correct  the  measured  line-intensity  ratio  for  self¬ 
absorption  effects.  Some  of  the  analytical  details  required  for  the 

(4  5) 

above  analysis  have  been  reported  by  Bader  ’  but  are  repeated  here  for 
the  sake  of  completeness. 

A  convenient  starting  point  is  the  solution  to  the  general 
equation  of  radiative  transfer  of  an  isothermal  slab  of  gas  of  thickness 
I  in  local  thermodynamic  equilibrium  in  which  scattering  effects  and  time- 
dependent  effects  have  been  neglected.  This  solution  takes  the  form^^*^^ 

=  B^(l  -  exp(-K^Jl))  (60) 

where 

I  =  Observable  spectral  intensity  at  frequency 

V  in  units  (energy/(tlme’area'solid  angle* 
frequency)) 

B  =  Blackbody  spectral  Intensity  at  frequency 

V  in  units  (energy/ (time 'area 'solid  angle* 
frequency)) 

=  Spectral  absorption  coefficient  including 
effects  of  induced  emission  in  units 
(1/length) 

I  =  Physical  thickness  of  gas  slab  in  units 
of  (length). 


The  blackbody  function,  B^,  takes  the  usual  form  of 

R  =  2hv^/c^ 

V  exp(hv/kT)  -1 


(61) 


where 


c  »  Speed  of  light 
h  ■  Planck  constant 
k  >  Boltzmann  constant 
T  ■  Temperature. 


The  Integrated  intensity  of  an  isolated  spectral  line  is 
obtained  by  integrating  Equation  (60)  over  frequency  v,  yielding 


^line 


/  I  dv  =  /  ■  exp(-K'i))  dv 

line  line 


(62) 


Planck's  function,  can  be  assumed  nearly  constant  for  the 
small  frequency  extent  of  the  spectral  line.  This  simplifies  Equation  (62) 
to  yield 

^line  °  ®v  dv  (63) 

o  line 


where  B^^  is  obtained  by  substituting  the  line  center  frequency,  v^,  into 
Equation  (61) . 

In  order  to  evaluate  the  Integral  in  Equation  (63) ,  a  form  for 

the  absorption  coefficient,  K^,  must  be  specified.  is  dependent  on 

the  broadening  mechanism  given  in  Section  III.  A  convenient  combination 

of  Doppler  and  collision-broadening  effects  (i.e.,  the  dominant  broadening 

mechanisms  in  the  RENT  environment)  can  be  modeled  into  a  composite 

absorption  coefficient  profile  known  as  a  Voight  profile.  For  this  profile, 

(21) 

the  analytical  form  of  is  given  by  Penner'  '  as 

K'  -  K'  H(a,5) 

V  o  ’ 

where 


H(a,C) 


a 

ir 


”  ,2. 

r  ■  ,exp(-y  ) _ 

J  2  2 

-<»  a  +  (C  -  y) 


dy 


K' 

o 

a 


/¥ 


1/2 

/SS’ 


/  K'dv 
line 


(65) 


(66) 

(67) 
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and 


■<  ' 


5  “ 


(V  -  v^) 

1/2  Y. 


/toT 


(68) 


Expressions  for  the  collision  and  Doppler  line  full  half-width  (i.e., 

Yp  and  y^)  are  given  in  Equations  (16)  and  (59),  respectively. 

The  Integrated  absorption  coefficient  given  in  Equation  (66) 
can  be  related  to  the  atomic  transition  probability  of  the  line  as' 

I  ,  ■  tA-  V  if  ^<1  -  )  .  (69) 

line  8nv  ®L 


Using  a  Boltzmann  distribution  to  describe  the  population,  N^,  of  the 
electronic  energy  levels  yields 


where 


St  8 

Nl  =  \  —  exp  (Eu  -  EL)/kT  =  ^  exp(hv^/kT) 

®u  °u 


(70) 


g^,  gj^  =  Degeneracy  of  upper  and  lower  states  of 
the  transition 


E^,  Ej^  =  Energy  of  upper  and  lower  energy  levels 
of  the  transition  (Joules) 

N^,  ~  Population  of  the  upper  and  lower  energy 

levels  (cm~^) 

=  Atomic  transition  probability  (sec  ^) 


and  the  copper  atom  number  density  in  the  upper  and  lower  states  can  be 
related  to  the  total  number  density  of  copper  atoms,  ,  using  the 
following: 


N  -  N  g  ^xp(-E^/kT) 

“  °Cu  P  exp(-E  /kT) 
i  ^  ^ 

exp(-E^ /kT) 

L  o^^  r  g^  exp(-E^/kT) 


(71) 


(72) 
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Hence,  substituting  Equation  (70)  and  (71)  into  Equation  (69),  the  integrated 
absorption  coefficient  can  be  written  as 


I 

line 


K'dv 

V 


8.  A 


j,  2  ®u  "uL  o 
OTTV  Cu 


-E  /kT  .  _ 

e  u  ,  hv  /kT 

- T -  (e  o  -1) 


o 

where  the  atomic  partition  function,  Q. ,  is  given  by 

A 

=  E  g^  exp(-E^/kT) 


(73) 


(74) 


The  so-called  "Voight  function"  (l.e.,  H(a,C))  in  Equation  (65)  is  not 

(22) 

solvable  in  closed  form.  However,  Armstrong  wrote  a  FORTRAN  IV  computer 
code  to  calculate  H(o,C)  to  an  accuracy  of  two  digits  in  the  sixth  signifi¬ 
cant  figure.  The  code,  as  reported  in  this  reference,  had  several  errors 
which  when  corrected  provided  a  convenient  means  of  computing  the  Voight 
function  for  use  in  determining  the  copper-line-intensity  profiles  of 
interest. 

Normalization  of  the  integrated  line  intensity  in  Equation  (11) 
by  Yjj  By^/  /in2  '  yields  the  following: 

I,  /ZnT  , 

-7— B -  =7/  (l-exp(-K'(l  H(a,C))  dC  .  (75) 

V,  line 

o 


The  integral  in  Equation  (75)  has  been  numerically  computed  for  various 

(23  ) 

values  of  and  o.  Janson  and  Korb'  give  this  parameter  in  tabular 
form  to  an  accuracy  of  four  significant  figures.  Use  of  this  tabulation 
results  in  the  determination  of  Integrated  line  intensities  to  an  accuracy 
of  0.40  percent  by  linear  interpolation  between  tabulated  values.  The 
table  covers  a  range  from  weak  to  fairly  strong  absorption. 

The  above  set  of  equations  and  numerical  procedures  represent 
the  line-intensity  profile  shape  and  the  Integrated  line  intensity  as  a 
function  of  pressure,  temperature,  copper  density,  and  optically-observed 
path  length  [l.e.,  requirements  for  Equation  (4)].  The  width  of  this 
emission  profile  can  also  be  determined  using  the  above  results.  Presently, 
no  closed-form  solution  is  available  to  represent  the  line  width,  y,  as  a 
function  of  the  above  dependent  variables  [l.e..  Equation  (5)  is  not 
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available  In  explicit  form].  However,  y(p,  T,  N^,  1)  is  available 
implicitly  from  the  above  computational  procedures.  The  formulations 
presented  in  this  section  can  therefore  by  used,  as  discussed  in  Section  I, 
to  correct  the  measured  line-intensity  ratio  for  self-absorption  effects. 

1.  Optically-Thin  Limit 

For  an  optlcally-thin  gas,  is  much  less  than  one  and  the 
exponential  term  in  Equation  (60)  can  be  approximated  by  the  first  two 
terms  in  the  series  expansion  for  the  exponential,  or 


exp(-K'Jl)  =  1  -  K'i 

V  V 


Substituting  Equation  (76)  into  Equation  (60)  yields 


I  ,  =  B  K'i 
V  V  V 


With  this  simplification.  Equation  (63)  can  also  be  rewritten  as 


line 


=  /  K’dv 

line 


(76) 


(77) 


(78) 


Finally,  substituting  Equations  (73)  and  (61)  into  Equation  (78),  we  find  that 

(79) 


^.T  exp(-E  /kT) 

'line  =  "'■o »o  ‘  - 

Cu 


For  the  optically— thin  case.  Equation  (79)  shows  that  the  total  intensity 
of  any  copper  emission  line  is  linearly  proportional  to  the  product  of 
copper  density  and  optical  path  length.  Therefore,  in  the  ratio  of  two 
optlcally-thin  lines,  these  factors  cancel  out  to  produce  the  following 
ratio  equation  using  v  «  1/X  to  obtain 


X2  ®Uj^  ^uLj^ 


I~  8  A  , 

2  1  ®U2  UL2 


exp 


rE  E 

-t  “1  -  '^27 


kT 


(80) 


where  and  X2  are  wavelengths  of  the  5106  and  5153  X  lines,  respectively. 
Solving  Equation  (80)  for  temperature  yields  the  desired  result 
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Using  the  above  gA  ratios,  plots  of  temperature  as  a  function  of  measured 
line  Intensity  ratio  for  each  of  the  above  cases  Is  shown  In  Figure  5. 

Note  the  large  dependence  of  the  ratio  small  changes  In  tempera¬ 
ture  for  values  less  than  AOOO  K.  This  curve  Indicates  that  the  region 

for  best  application  of  the  5106/5153  Intensity  ratio  Is  between  4000  and 
8000  K.  Also  note  the  Importance  of  having  an  accurate  value  for  the 

transition  probability  ratio  of  a  given  line  pair.  A  7  percent  error 

would  have  resulted  In  the  measured  temperature  under  optlcally-thln 
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TABLE  4 


conditions  if  Corliss  and  Bozman's  gA  values  were  utilized  near  6000  K 
Instead  of  Calspan's  results. 


line  Intensities  is  a  function  only  of  the  gas  temperature  (i.e.,  copper 
density  and  path  length  cancel  out).  Calculations  of  the  optically-thick 
line-intensity  ratio  were  performed  using  equations  and  constants  given  in 
Sections  III  and  IV.  The  results  are  also  shown  plotted  in  Figure  5  along 
with  the  optically-thin  computations.  Note  that  significant  temperature 
error  could  result  if  measurements  of  optically-thick  lines  are  used  in 
conjunction  with  optically-thin  theory.  The  opposite  is  also  true; 
therefore,  the  copper  density  must  be  known  or  large  errors  are  possible 
utilizing  the  line-ratio  method. 
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SECTION  V 

DESCRIPTION  OF  COMPUTATIONAL  PROCEDURES 

Analytical  foraulatlons  for  computing  the  line  intensity  pro¬ 
files  of  the  5106  and  5153  A  copper  lines  which  are  given  in  Sections  III 
and  IV  were  coded  into  a  FORTRAN  IV  digital  computer  program.  For  a 
detailed  description  of  this  program,  the  reader  should  consult  Reference  25. 

The  input  variables  required  by  the  code  are  as  follows: 

Gas  pressure,  p,  in  atm 
Gas  temperature,  T,  in  “K 
Optical  path  length,  8,,  in  cm 

3 

Copper  number  density,  in  atoms/cm 

Various  atomic  constants. 

The  code  first  calculates  the  true  line  shape  (i.e.,  inf initely-narrow 
Instrument  apparatus  function)  using  Equations  (60),  (61),  (64)  through 
(68),  and  (73).  The  results  of  this  calculation  are  plotted  along  with  a 
printout  of  the  ratio  of  line  peak  intensities,  the  integrated  line- intensity 
ratio,  and  the  ratio  of  the  widths  of  the  two  lines  (i.e.,  Yj^/Y2)’  ^  table 

look-up  routine  is  supplied  for  the  thermodynamic  calculations  to  provide 
values  for  the  number  density  of  neutrals  and  electrons  given  the  tempera¬ 
ture  and  pressure.  The  table  was  constructed  using  data  provided  in  Neel's 
computations^  .  Table  look-up  results  are  also  printed  out  as  part  of 
the  code  output.  The  partition  function,  Q^,  is  also  provided  as  a  table 
look-up  based  on  calculations  using  Equation  (74). 

Computation  of  true  line  shape  is  accomplished  by  calculating 

(22) 

the  "Volght  function"  using  Armstrong’s  numerical  methods  '  lor 
wavelength  positions  from  the  line  center  to  a  wavelength  where  the 
intensity  was  less  than  0.1  percent  of  the  peak  intensity  value.  This 
line  contour  is  then  numerically  analyzed  to  determine  the  peak  intensity 
and  the  full  half-width  of  the  line.  The  profile  was  also  numerically 
Integrated  using  Simpson's  rule  to  obtain  the  integrated  line  intensity 
value.  These  true  line-shape  results  are  printed  out  along  with  the  input 
parameters  for  the  particular  calculations. 

Since  most  spectral  data  are  obtained  with  a  finite  apparatus 
function  width,  some  convolution  occurs  and  the  true  line  shape  is  altered 
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due  to  the  shape  of  the  apparatus  function.  Therefore,  the  code  also 
Included  a  convoluting  routine  to  convert  the  above  true  line  contour  into 
the  profile  actually  measured  by  the  spectroscopic  equipment.  The 
apparatus  function  used  to  perform  this  convolution  was  assumed  to  be 
Gaussian  in  nature.  Hence,  the  function  was  completely  defined  for  a 
given  full  half-width  of  the  Gaussian  line  profile  which  is  the  final 
input  parameter.  This  apparatus  function  width,  y was  determined  from 
experimental  data  obtained  using  an  infinltely-thln  calibration  line  from 
standard  sources  (i.e.,  copper  hollow-cathode  lamp).  A  provision  was  also 
incorporated  into  the  code  to  allow  separate  apparatus  functions  for  each 
of  the  two  lines  of  interest  (i.e.,  and  Ya2^"  Other  apparatus  function 
shapes  could  be  coded  if  necessary.  Once  the  convolution  calculation  was 
performed,  the  convoluted  profile  was  also  numerically  analyzed  to  provide 
peak  intensity  ratios,  integrated  intensity  ratios,  and  apparent  line 
widths  for  actual  measurement  conditions.  These  results  were  also  printed 
as  part  of  the  code  output.  Measured  spectral  data  could  then  be  compared 
directly  with  theoretical  values  to  avoid  the  necessity  for  deconvolution 
of  the  measured  line-intensity  profile  data. 

Several  other  parameters  are  printed  out  for  each  line  of 
interest.  They  are  as  follows: 

K^,  a  -  From  Equations  (66)  and  (67),  respectively 
-  Equation  (61) 

-  Peak  line  intensity 

Iline  “  Integrated  line  intensity 

Y  -  Actual  full  half-width  of  the  line  in  X 

Yd»  ~  Doppler,  Stark,  and  Van  der  Waals  con¬ 

tributions  to  the  true  width,  in  X. 

The  above  code  output  completely  defines  the  shape  of  and  intensity  param¬ 
eters  for  the  5106  and  5153  X  lines.  Pressure,  optical  path  length,  and 
apparatus  function  widths  are  usually  held  constant  while  temperature  and 
copper  density  are  varied  to  generate  a  matrix  of  dependent  variables  using 
the  above  code.  Plots  of  these  dependent  parameters  as  a  function  of 


temperature  and  copper  density  can  be  used  to  reduce  spectroscopic  data 
from  the  arc-tunnel  flow  In  order  to  obtain  the  gas  flow  static  tempera¬ 
ture  and  unknown  copper  density  corrected  for  self-absorption  effects  on 
the  emission  spectra. 


SECTION  VI 

DISCUSSION  OF  ANALYTICAL  RESULTS 


The  computer  code  described  briefly  in  Section  V  was  exercised 
for  the  following  input  conditions: 

p  =  80  atm 

H  =  2  cm 

=  1.64  A 

Ya2  =  1.41  X 

T  =  4500-8000  K  in  500  K  increments 

Nq  =  10^^,  10^^,  2  X  10^^,  5  X  10^^,  10^^,  2  x  10^^, 

1 A  Y 

5  x  lO'*'”  atoms/cm  . 

It  should  be  noted  that  the  present  choice  of  apparatus  function 
widths  was  not  fortuitous,  but  represents  values  obtained  in  later  experi¬ 
mental  work  by  the  present  author.  The  two  functions  had  different 
half-widths  due  to  the  slight  difference  in  the  focal  properties  of  the 
measurement  instrument.  These  apparatus  function  widths  also  are  representa¬ 
tive  of  values  used  by  Bader^^^  in  a  four-channel  polychromater  instrument 
which  employed  a  trapezoidal  apparatus  function  that  had  a  full  half-width 
of  1.79  X  and  a  base  width  of  2.66  X.  Therefore,  the  observations  stated 
below  have  validity  for  both  past  and  present  experimental  efforts.  Results 
of  theoretical  computations  for  true  line  contours  (l.e.,  =  0)  and 

relative  line  intensities  are  shown  in  Figure  6  for  T  =  6000  K  as  a  func¬ 
tion  of  flow  copper  density.  In  Figure  7,  the  same  results  are  shown  for 
both  lines  but,  in  this  case,  the  true  line  contours  have  been  convoluted 
with  the  proper  full  half-width,  Gausslan-shaped  apparatus  functions  to 
indicate  measured  line  profile  conditions. 

Note  that  at  fixed  temperature  and  pressure  (i.e.,  constant 
colllslonal  and  Doppler  broadening) ,  the  following  effects  occur  for  each 
of  the  lines  as  the  copper  density  increases. 

•  Line  full  half-width  Increases  due  to  self¬ 
absorption  effects  near  the  center  of  each  line. 
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Wavelength  From  Line  Center,  A 


TRUE  LINE  CONTOUR  AND  RELATIVE  LINE  INTENSITY  CHANGES  WITH 
INCREASING  COPPER  NUMBER  DENSITY  (p  -  80  atm,  i  -  2  cm. 


Ratotive  Intensity 


•  The  width  of  the  5106  X  line  is  less  than  the 
5153  I  line  at  optically-thin  conditions  but  the 
width  of  the  former  line  increases  faster  with 
increasing  copper  density  until  the  two  lines  are 
nearly  of  equal  width  at  high  copper  density. 

•  Relative  line  intensity  changes  such  that  the 
5106  to  5153  X  peak  intensity  ratio  decreases 
until  the  lines  are  nearly  equal  in  peak  intensity 
at  high  copper  density. 

•  Integrated  intensity  ratios  follow  the  same  trend 
as  the  peak  ratios  but  to  a  lesser  degree. 

•  As  the  optical  thickness  Increases,  the  lines 
become  increasingly  flattened  (i.e.,  constant 
intensity)  near  the  line  center. 

The  effect  of  a  relatively  broad  apparatus  function  on  the 
spectral  line  contours  is  shown  in  Figure  7.  Note  that  details  of  self¬ 
absorption  effects  are  masked  somewhat  by  the  broad,  unequal  apparatus 
functions,  but  that  the  overall  changes  in  shape,  width,  and  peak 
intensity  of  the  lines  noted  above  can  still  be  recognized. 

In  Figures  6  and  7,  the  magnitude  of  self-absorption  broadening 
with  increased  copper  density  can  be  readily  observed  in  comparison  to 
the  ordinary  broadening  mechanisms  which  dominate  at  low  copper  density. 

Note  that  at  a  copper  density  of  approximately  2  x  10^^  atoms/cm^,  the 
true  line  width  of  the  5106  A  line  due  to  self-absorption  effects  is  about 
twice  as  wide  as  it  would  be  for  optically-thin  conditions. 

Examination  of  the  spectra  in  Figures  6  and  7  also  brings  into 
sharp  focus  the  fact  that  large  errors  could  be  encountered  if  a  simple 
ratio  of  Integrated  or  peak  intensity  ratios  was  used  to  infer  gas  tempera¬ 
ture  without  a  knowledge  of  the  flow  copper  density.  However,  examination 
of  the  plots  in  these  figures  also  Indicates  that  sufficient  changes  occur  in 
the  line  profile  shape  and  full  half-width  as  copper  density  Increases  to 
imply  that  both  copper  density  and  temperature  could  be  simultaneously 
determined  from  a  measurement  of  the  entire  line  intensity  profile  (i.e., 
peak  intensity  and  line  full  half-width) . 


To  accomplish  this  simultaneous  temperature/copper  density 
measurement.  It  Is  necessary  to  first  compute  line  peak  Intensity  ratio 
and  line  full  half-width,  for  a  series  of  possible  copper  density/temperature 
combinations.  Results  of  these  computations  are  shown  In  Figures  8  and  9 
for  an  Inf Inltely-thln  apparatus  function  as  well  as  one  which  represents  a 
more  realistic  set  of  Instrument  full  half-widths  (l.e.,  Ya^=1>64  A, 

YA2  "  1-^1  A)- 

The  particular  values  of  plotting  parameters  shown  In  Figures  8 
and  9  were  chosen  for  several  reasons.  Peak  Intensity  ratios  were  utilized 
rather  than  Integrated  ratio  values  In  order  to  simplify  the  data  reduction 
procedure.  Absolute  magnitude  of  the  5106  A  line  width  was  utilized  rather 
than  a  ratio  of  the  two  line  widths  since  Calspan  data  determined  the 
broadening  constants  for  the  5106  A  line  only. 

The  curves  In  Figure  9  could  be  utilized  directly  to  reduce 
spectroscopic  data  taken  with  the  same  apparatus  function.  However,  a 
cumbersome  Iteration  process  would  be  required  whereby  measured  values 
of  Intensity  ratio  and  5106  A  line  width  are  used  to  locate  a  set  of 
temperature  and  copper  density  values  that  are  consistent  with  both  of 
the  measured  parameters.  A  simpler  way  to  reduce  the  data  can  be  obtained 
if  cross-plots  are  constructed  from  the  theory  plots  in  Figure  9.  These 
cross-plots  are  obtained  by  picking  constant  values  of  intensity  ratio  and 
line  width  and  finding  the  points  of  intersection  of  these  values  with  the 
constant  temperature  and  copper  density  curves  in  Figure  9.  New  curves 
can  then  be  constructed  with  temperature  and  density  for  the  variables  as 
shown  in  Figure  10. 

Again  in  Figure  10,  the  importance  of  knowing  the  copper  density 
is  evident  in  order  to  infer  temperature  from  line-intensity  ratio  data. 

Note  that  only  for  narrow  (l.e.,  optically-thin)  lines  is  the  line  ratio 
Independent  of  flow  copper  density. 

Comparison  of  Figures  8  and  9  also  Indicates  that  a  greater 
accuracy  would  be  obtainable  for  this  data  reduction  process  if  the  plots 
In  Figure  10  were  constructed  with  calculations  using  a  narrow  apparatus 
function  (l.e.,  Ya  <<  A).  This  is  especially  true  for  the  line  width 

A 

data  because  at  moderate  copper  densities  the  spread  between  the  constant 
temperature  lines  is  significantly  greater  for  “  0  than  for  -  1.5  A. 
However,  for  high  copper  densities  (l.e.,  broad  true  line  widths),  even 
the  broader  apparatus  function  is  quite  acceptable. 
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5106  A  Line  Full  Half- Width, 


Temperoture,  “K 


5106  A  Line  Full  Half -Width, 


5106/5153  Peak  Intensity  Ratio 


Temperature ,  ®K 


In  Figure  10,  the  intersection  of  the  two  families  of  curves 
defines  compatible  temperature/copper  density  combinations  (l.e.,  solutions 
to  the  desired  data  reduction  process).  Hence,  Figure  10  can  be  used  to 
quickly  determine  which  temperature  and  copper  density  combination  that  best 
agrees  theoretically  with  the  measured  parameters.  The  process  of  reducing 
the  data  Is  as  follows: 

(1)  The  curve  which  most  closely  represents  the 
measured  5106  X  line  width,  Yj^»  is  first  located 
in  Figure  10  using  linear  interpolation,  if 
necessary. 

(2)  The  corresponding  curve  representing  the 
measured  peak  intensity  ratio  (Ij^/12)  is  then 
located  as  above. 

(3)  The  intersection  of  these  two  curves  gives  the 
required  copper  density  and  temperature. 

A  more  accurate  method  of  analysis  would  involve  construction 
of  curves  such  as  those  in  Figure  10  for  many  more  constant  width  and 
Intensity  ratio  values.  Note  also  that  the  above  data  reduction  procedure 
could  be  easily  mechanized  and  coded  into  a  computer  program  to  automate 
the  data  reduction  process.  This  automation  was  accomplished  in  the  present 
study  and  the  code  was  available  for  computer  data  reduction. 


SECTION  VII 
ERROR  ESTIMATE 


The  text  of  this  report  contains  a  method  to  correct  the  simple 
line-intensity  ratio  temperature-measurement  technique  for  self-absorption 
effects.  However,  a  discussion  of  this  correction  method  would  not  be 
complete  without  some  indication  as  to  the  accuracy  of  the  correction 
procedure.  First  of  all,  it  is  necessary  to  point  out  that  several  other 
assumptions  must  be  valid  before  the  line-ratio  method  will  provide  the 
most  accurate  results  possible  even  for  optically-thin  conditions.  These 
assumptions  are 

•  The  temperature,  pressure,  and  copper  density  are 
assumed  to  be  spatially  uniform  over  the  emission 
region;  that  is,  the  above  properties  must  be 
spatially  homogeneous  over  the  optical  path  length. 

•  A  corollary  to  the  above  assumption  is  required 
when  radiation  from  a  model  nose  tip  is  being  ex¬ 
amined  to  infer  model  stagnation  temperature. 

This  corollary  states  that  the  free-stream  must 
emit  and  absorb  a  negligible  amount  of  radiation 
compared  to  these  effects  in  the  nose-tip  region. 

Further,  it  is  assumed  that 

•  Local  thermodynamic  equilibrium  (LTE)  exists  such 

that  population  of  the  electronic  mode  of  copper 

is  represented  by  a  Boltzmann  distribution  described 

d 

by  some  temperature,  T  ,  that  is  equal  to  the 

el 

static  temperature  of  the  gas  (i.e.,  T^,^  =  T) . 

•  The  final  assumption  is  that  the  measurement  time 
is  short  compared  with  any  temporal  fluctuations 
in  flow  properties. 

We  will  presume  that  the  above  assumptions  are  correct  such 
that  the  only  error  source  is  considered  to  be  due  to  self-absorption 
effects.  The  question  then  Is,  "How  accurate  Is  the  self-absorption 
correction  procedure?" 
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Several  dominant  sources  can  be  listed  as  possible  causes  of 
error  In  the  self-absorption  analysis  procedure,  and  they  can  be  grouped 
Into  three  categories: 

•  Measurement  Inaccuracies 

•  Erroneous  thermodynamic  assumptions 

•  Faulty  analytical  assumptions. 

Measurement  Inaccuracies  are  always  possible  and  Involve  the  Inability  to 
determine  the  exact  value  of  the  experimental  line  width  and  Intensity 
ratio.  Without  valid  measurement  values  In  the  first  place,  It  Is  Impossible 
to  accurately  correct  for  self-absorption  errors.  Another  measurement 
error  source  Is  the  determination  of  an  exact  apparatus  function  and  Its 
associated  Gaussian  full  half-width.  Also  an  accurate  value  for  the  optical 
path  length,  I,  must  be  obtained  from  photographic  observation  of  the 
radiating  gas.  Finally,  the  accuracy  of  measurements  made  by  other  experi¬ 
menters  can  lead  to  errors  In  the  present  analysis  since  atomic  constants 

(l.e.,  g, A, /g-A„,  C,'s)  determined  by  other  authors  are  used  In  the  present 
1  1  Z  Z  D 

correction  procedure. 

Regarding  the  thermodynamic  assumptions  used  In  the  data  re¬ 
duction  procedure.  It  should  be  noted  that  the  gas  pressure  Is  assumed  to 
be  known  from  flow  measurements  obtained  by  AFFDL  personnel.  The  pressure 
Is  also  assumed  to  be  constant  over  the  entire  wind-tunnel  run.  Further, 

It  Is  assumed  that  the  gas  Is  In  chemical  equilibrium  such  that  the  number 
density  of  neutrals  and  electrons,  which  are  used  In  the  line-broadening 
calculations,  can  be  determined  from  the  known  pressure  and  the  assumed 
temperature. 

The  final  error  category  Is  the  Inaccuracy  of  theoretical  models 
used  to  predict  the  line  shape  and  Intensity.  It  Is  felt  that  the  "Volght" 
line  profile  Is  quite  accurate  and  that  Doppler  broadening  Is  well  under¬ 
stood.  On  the  other  hand,  the  Llndholm  theory  for  collision  broadening  may 
not  be  a  true  representation  of  broadening  at  the  high  pressure  encountered 
In  the  RENT  arc-heated  flow.  However,  It  Is  difficult  to  assign  an  error 
value  to  this  type  of  possible  Inaccuracy  because  there  Is  presently  no 
exact  theory  to  compare  with  the  Llndholm  computational  results.  Perhaps 
the  line  shift  Information  obtained  in  the  RENT  environment  will  sufflcently 


corroborate  the  application  of  Llndholm's  model  to  the  arc-heated  flow 
environment.  The  theory  as  coded  does  predict  the  width  and  intensity 
ratio  that  were  measured  by  Calspan  in  shock-tube  experiments.  These 
latter  measurements  were  made  at  conditions  similar  to  the  E5ENT  environ¬ 
ment  so,  in  actuality,  the  subtraction  of  the  optically-thin  line  width 
from  the  measured  line  width  to  infer  copper  density  tends  to  be  empirically 
correct . 

The  accuracy  of  this  self-absorption  correction  procedure  was 
determined  by  varying  the  parameters  whose  values  cannot  be  specified  with 
certainty  by  a  fixed  amount  about  some  base  value.  Plots  such  as  those 
in  Figure  10  could  have  been  constructed  for  each  parameter  variation  and 
the  corresponding  variation  in  corrected  temperature  could  be  noted.  This 
tedious  but  straightforward  procedure  provides  an  indication  of  the 
sensitivity  of  the  inferred  temperature  to  inaccuracies  in  the  above  data 
reduction  parameters.  Then,  an  estimate  could  be  made  of  the  possible 
inaccuracy  in  these  parameters  for  the  actual  conditions  encountered  in  a 
given  experiment.  Finally,  the  results  from  the  sensitivity  data  could 
then  be  coupled  to  the  above  error  estimates  to  compute  a  possible  error 
in  the  corrected  temperature  based  on  all  the  possible  inaccuracies  in  the 
data  reduction  parameters. 

Actually,  this  process  should  have  been  repeated  for  each  experi¬ 
mental  data  point  because  the  resulting  temperature  correction  error  would 
vary  with  the  particular  region  of  Figure  10  where  the  data  fell.  However, 
a  representative  error  was  obtained  by  examining  a  data  point  that  falls 
nearest  the  average  temperature  and  copper  density  values.  The  baseline 
intensity  ratio  of  1.65  and  the  5106  line  width  of  2.69  A  were  taken  from 
actual  arc-tunnel  data  that  represented  a  temperature  and  copper  density 
near  the  average  for  all  the  data  obtained  (i.e.,  ~  ^^38  K  and 

)u  n  ”  1.10  X  10^^  atoms/cm^). 

Cu  baseline  ' 

The  above  error  estimate  process  was  then  undertaken  using 
parameters  employed  to  produce  Figure  10  as  the  baseline  case.  Each 
Important  parameter  was  then  varied  about  this  baseline  by  a  prescribed 
amount  and  the  corresponding  temperature  variation  was  calculated.  Some 
of  the  analytical  parameters  have  a  small  effect  on  the  computations  so 
they  were  left  out  of  the  analysis.  The  computer  program  written  to 
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automatically  perform  the  required  Iterations  was  utilized  such  that  no 
Interpolation  or  replottlng  process  was  necessary.  This  scheme  provided 
more  accurate  sensitivity  estimates. 

Results  of  this  error  analysis  are  shown  In  Table  5.  For  each 
selected  variation,  APar,  In  the  baseline  parameter,  a  corresponding  Inferred 
temperature  variation,  AT,  was  computed.  Therefore,  a  sensitivity  parameter, 
S,  can  be  computed  where 


_ B  _  Error  In  Temperature 

APar/Par„  Error  In  Parameter 


This  parameter,  S,  represents  the  sensitivity  of  the  Inferred  temperature 

to  errors  or  uncertainties  In  the  value  of  a  given  parameter.  Hence,  large 

sensitivities  coupled  with  large  parameter  uncertainties  result  In  large 

measurement  errors.  The  results  In  Table  5  Indicate  that  the  sensitivity 

Is  small  for  all  the  parameters  of  interest. 

Also  given  in  Table  5  is  an  estimate  of  the  actual  uncertainty 

in  a  given  parameter,  e  .  These  estimates  represent  conservative  values 

par 

from  the  author's  own  data  as  well  as  that  of  other  authors.  A  value  of 
50  percent  accuracy  of  for  the  5153  line  was  based  on  the  fact  that  no 
line  width  measurements  were  made  for  this  line  but  the  measurements  for 
the  other  line  were  within  this  amount  of  error  from  the  theoretical  values. 

Using  these  error  estimates,  the  actual  temperature  error,  e^, 
for  each  of  these  error  sources  could  be  computed  using 

Temperature  Error  =  e„  =  S  •  e  .  (91) 

T  par 

These  errors  are  also  shown  in  Table  5  and  most  are  less  than  1  percent. 

The  largest  error  source  is  due  to  the  uncertainty  in  the  ratio  of  gA  factors 
for  the  two  lines.  Uncertainty  in  Cg's  is  also  important.  The  total 
possible  error  equals  13.5  percent  (i.e.,  ±7  percent),  which  is  a  simple 
sum  >f  all  the  error  sources  in  label  5.  Therefore,  we  can  write  that  the 
measured  stagnation  temperature,  T^,  is  given  by 

T^  =  6538  ±  458  K  .  (92) 

The  resultant  inferred  enthalpy  error  is  larger  because  the  sensitivity 
of  enthalpy  change  to  temperature  change  is  greater  than  one.  Hence,  the 
measured  stagnation  enthalpy  was 
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Atoms/cm 


H  -  5400  ±  625  [Btu/lb  ] 


or  the  error  is  ±12  percent  for  self-absorption  effects  only.  This 
enthalpy  measurement  error  could  be  reduced  to  less  than  ±10  percent  by 
judicious  reduction  of  certain  error  sources  listed  in  Table  5,  especia 
those  relating  to  uncertainties  in  the  atomic  constants. 


SECTION  VIII 
SUMMARY  AND  CONCLUSIONS 

The  text  of  this  report  contains  detailed  analytical  formulae 
and  appropriate  empirical  constants  required  to  predict  shape  and  intensity 
of  the  copper  emission  spectrum  from  a  high-temperature,  high-pressure 
plasma  including  self-absorption  effects.  In  particular,  the  spectrum  of 
the  5106  and  5153  A  copper  lines  was  calculated  for  conditions  typical  of 
the  gas  flow  environment  in  the  Reentry  Nose-Tip  (RENT)  Facility  at  the 
Air  Force  Flight  Dynamics  Laboratory  (AFFDL) .  In  previous  work  at  AFFDL, 
a  simple  intensity  ratio  from  the  above  lines  was  used  to  infer  gas 
temperature  in  the  RENT  facility.  However,  the  arc-heated  flow  contains 
an  unknown  amount  of  copper  contaminant,  so  the  above  theoretical  formulation 
was  required  in  order  to  provide  a  method  of  correcting  the  simple  line 
ratio  measurement  for  self-absorption  effects. 

Development  of  the  analytical  expressions  presented  in  this 
report  includes  a  discussion  of  various  line-broadening  mechanisms  and 
an  estimate  of  the  relative  magnitude  of  each  of  the  various  operative 
broadening  effects.  The  empirical  constants  required  to  predict  collisional 
line  broadening  of  the  5106  A  line  for  optically-thin  conditions  were 
derived  from  spectral  data  measured  by  Calspan  under  a  contract  to  AFFDL. 
Data  required  to  compute  the  magnitude  of  other  broadening  effects  were  not 
provided  by  this  source  and  had  to  be  obtained  from  other  references  in 
the  open  literature. 

Based  on  these  calculations,  it  was  found  that  the  dominant 
broadening  mechanisms  in  the  RENT  environment  included 

•  Self-absorption  broadening  due  to  excessive 
amounts  of  copper  contaminant 

•  Van  der  Waals  broadening  due  to  collisions 
between  emitting  atoms  and  neutral  particles. 

Other  less  important  broadening  phenomena  that  were  included  for  complete¬ 
ness  were  Doppler  and  Stark  broadening  effects.  All  other  broadening 

mechanisms  were  found  to  be  negligible  for  RENT  flow  conditions  (i.e,, 

16  3 

pressure  *  100  atm,  temperature  =  6000  K,  copper  density  =  10  atoms/cm  ) . 
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Dominance  of  colllslonal  broadening  mechanisms  leads  to  the 
use  of  a  "Voight"  Line  contour  shape.  This  line  profile  was  used  to 
predict  the  line  width  and  peak  line-intensity  ratio,  as  well  as,  inte¬ 
grated  line-intensity  ratio  as  a  function  of  flow  copper  density,  gas 
pressure,  optical  path  length,  and  gas  temperature.  Standard  line  intensity 
theory  described  by  Penner  was  used  in  these  computations  with  transition 
probability  data  also  determined  by  Calspan. 

Results  of  these  line  shape  and  intensity  ratio  calculations 
showed  that,  for  conditions  typical  of  the  RENT  flow  environment,  self¬ 
absorption  of  copper  line  radiation  was  significant  such  that  a  large  error 
in  measured  temperature  would  be  encountered  if  the  simple  line-intensity 
ratio  was  not  corrected  for  self-absorption  effects.  It  was  also  shown 
that  correction  of  the  simple  line-ratio  method  was  indeed  possible  if  a 
simultaneous  measurement  was  made  of  the  full  half-width  of  the  5106  A 
copper  line.  These  two  measurements  were  shown  to  provide  the  tedious, 
but  necessary,  means  to  correct  the  simple  line-ratio  method  such  that  a 
more  accurate  temperature  could  be  inferred  from  spectroscopic  emission 
data.  The  present  estimate  of  the  temperature  error  resulting  from  un¬ 
certainties  in  this  self-absorption  correction  procedure  is  ±7  percent  for 
typical  RENT  peaked-profile  stagnation  flow  conditions.  This  corresponds 
to  a  ±12  percent  error  in  the  enthalpy  inferred  from  the  self-absorbed 
copper  line  spectral  data.  This  self-absorption  correction  error  could 
be  reduced  if  more  accurate  atomic  constants  were  available.  Other  error 
sources  also  need  to  be  examined  before  a  final  error  estimate  can  be 
assigned  to  the  copper-llne-ratlo  techniques  for  measuring  enthalpy  in  an 
arc-heated  flow. 
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